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Abstract. We establish the integrability of the last open case in the Kozlov- 
Treshchev classification of Birkhoff integrable Hamiltonian systems. The tech- 
nique used is a modification of the so called quadratic Lax pair for D n Toda 
lattice combined with a method used by M. Ranada in proving the integrability 
of the Sklyanin case. 



1. Introduction 

In this paper we prove the complete integrability of the Hamiltonian 
system defined by 

n 1 n—1 

(1) H = J2 ~P 2 i + Y1 e " l ~ qi+1 + e qn ~ 1+qn + e~ qi + e- 2qi . 

i=l i=l 

The integrability of this system was conjectured in [16] and in the 
book of V.V. Kozlov [17]. This system appears first in the classification 
of Birkhoff integrable systems by Kozlov and Treshchev [16]. The clas- 
sification involves systems with exponential interraction with sufficient 
number of integrals, polynomial in the momenta. The classification 
gives necessary conditions for a system with exponential interraction 
to be Birkhoff integrable. The integrability (or not) of each system 
in the list was established case by case using various techniques. The 
only open case which remains is the case of system ([T]). For this rea- 
son, this last case became sort of famous and we refer to it as the last 
integrable case of Kozlov- Treshchev potentials. We now give a brief his- 
torical review of this area, including previous progress in establishing 
the integrability of Birkhoff integrable systems. 

We begin with the following general definition which involves systems 
with exponential interaction: Consider a Hamiltonian of the form 

1 N 

(2) H = -(p,p) + 5>^ , 

i=i 

where q = (qx, . . . , q n ), p = (pi, . . . ,p n ), Vi, ■ ■ ■ , v N are vectors in R n 
and ( , ) is the standard inner product in R n . Following [17] we call 
the set of vectors A = {vi, . . . , Vjv} the spectrum of the system. 
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Let M be the N x N matrix whose elements are 

Mij = (vi,Vj) . 

Hamilton's equations of motion can be transformed by a generalized 
Flaschka transformation to a polynomial system of 2N differential 
equations. The transformation is defined as follows: 

(3) ai = -e^\ 6 i = (v i ,p). 

We end-up with a system of polynomial differential equations 

Ofc = akbk 

N 

(4) b k = M ^ ai ■ 

1=1 

Equations (jlj) admit the following two integrals 

N N 

(5) F 1 = ]TAA, F 2 = l[a^ 

i=l i=l 

provided that there exist constants Aj such that X/i=i -^ v « = 0- Such 
integrals always exist for N > n. One can define a canonical bracket 
on the space of variables (a*, &*) by the formula 

= (Vi,v 3 -)ai 

and all other brackets equal to zero. The integrals F% and F 2 are 
Casimirs of this bracket. 

An interesting special case of (j2J) occurs when the spectrum is a 
system of simple roots for a simple Lie algebra Q . In this case N = I = 
rankC?. It is worth mentioning that the case where N, n are arbitrary 
is an open and unexplored area of research. The main exception is 
the work of Kozlov and Treshchev [16J where a classification of system 
(T5]) is performed under the assumption that the system possesses n 
polynomial (in the momenta) integrals. We also note the papers by 
Ranada [22J, Annamalai, Tamizhmani [3], Emelyanov [9], Emelyanov 
and Tsygvintsev pi)]. Such systems are called Birkhoff integrable. For 
each Hamiltonian in (j2J) we associate a Dynkin type diagram as follows: 
It is a graph whose vertices correspond to the elements of A. Each pair 
of vertices Vj, Vj are connected by 

4(v„v J ) 2 

(VijViXyj, Vj) 

edges. 
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Example 1. The origin of systems of exponential interraction is the 
classical Toda lattice which corresponds to a Lie algebra of type A n _i. 
In other words N = I = n — 1 and we choose A to be the set: 

Vl = (1,-1,0,...,0), ,v„_i = (0,0,..., 0,1,-1) • 

The graph is the usual Dynkin diagram of a Lie algebra of type A n _i. 
The Hamiltonian becomes: 

n j n— 1 

(6) H(q 1} ...,q n ,Pi,...,Pn)=J2 2 p2i+ ^2 e " +1 ' 

i=l i=l 

which is the well-known classical, non-periodic Toda lattice. The in- 
terability of this system was established in [11], [12], [H], [18], [T9] . 

As we already mentioned, the Toda lattice was generalized to the case 
where the spectrum corresponds to a root space of an arbitrary sim- 
ple Lie group. This generalization is due to Bogoyavlensky [I]. These 
systems were studied extensively in [15] where the solution of the sys- 
tems was connected intimately with the representation theory of simple 
Lie groups. There are also studies by Olshanetsky and Perelomov [20] 
and Adler, van Moerbeke [2]. The case of D n -Toda lattice plays an 
important role in the present paper. 

It is more convenient to work, instead with the space of the natural 
(q,p) variables, with the Flaschka variables (a, b) which are defined by: 



(?) a, = ie^) i = l,2,...,N 

h = ~\Vi i = 1,2, . . . ,n . 

We end-up with a new set of polynomial equations in the variables 
(a, b). One can write the equations in Lax pair form, see for example 
[2T] . The Lax pair (L(t), B(t)) in Q can be described in terms of the 
root system as follows: 

l i 

i=l i=l 
I 

B(t) = ^a,(t)(e ai - e_ Q J . 
i=i 

As usual h ai is an element of a fixed Cartan subalgebra and e ai is a root 
vector corresponding to the simple root «j. The Chevalley invariants 
of Q provide for the constants of motion. In this paper we make use of 
transformation ([7]) as opposed to transformation ([3]). 

The first important result in the search for integrable cases of system 
(T5]) is due to Adler and van Moerbeke pp. They considered the special 
case where the number of elements in the spectrum A is n + 1 (i.e., 
iV = n+ 1). Furthermore, they made the assumption that any n vectors 
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in the spectrum are independent. Under these conditions a criterion 
for algebraic integrability is that 



should be in the set {0,-1,-2,...} for all i ^ j. The method of 
proof in [TJ is based on the classical method of Kovalevskaya. The 
classification obtained corresponds to the simple roots of graded Kac- 
Moody algebras. The associated systems are the periodic Toda lattices 
of Bogoyavlensky [4 J . The complete integrability of these systems using 
Lax pairs with a spectral parameter was already established in [2]. 

The next development in the study of system (j2J) is the work of 
Kozlov and Treshchev on Birkhoff integrable systems. A system of the 
form (j2J) is called Birkhoff integrable if it has n integrals, polynomial 
in the momenta with coefficients of the form 

^A ie <u ^ q> , A 3 G R, u 3 GR", 

whose leading homogeneous forms are almost everywhere independent. 
We remark that in the definition given in the book of Kozlov [T7] there 
is no assumption on involutivity of the integrals. In [16j it is proved 
that the polynomial integrals are in involution. For this reason, we will 
not deal with the involution of the integrals; it follows from Lemma 5, 
p. 567 of [16]. The terminology has its origin in the work of Birkhoff 
who studied the conditions for the existence of linear and quadratic 
integrals of general Hamiltonians in two degrees of freedom. A vector 
in A is called maximal if it has the greatest possible length among all 
the vectors in the spectrum having the same direction. Kozlov and 
Treshchev proved the following theorem: 

Theorem 1. Assume that the Hamiltonian (TJ|) is Birkhoff integrable. 
Let Vj be a maximal vector in A and assume that the vector Vj G A is 
linearly independent o/vj. Then 

2(vi, Vj) 

(Vi, Vi) 

lies in the set {0, —1, —2, . . . }. 

Note that the condition of the theorem is exactly the same as con- 
dition (jSJ) of Adler and van Moerbeke. Of course theorem 1 is more 
general since there is no restriction on the integer iV (the number of 
summands in the potential of (j2J)). It turns out, however, that N can- 
not be much bigger than n. In fact, it follows from the classification 
that N < n + 3. A system of the form (j2J) is called complete if there 
exist no vector v such that the set A U {v} satisfies the assumptions 
of Theorem [TJ In [16] there is a complete classification of all possible 
Birkhoff integrable systems based on Theorem [TJ The Dynkin type 
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diagram of a complete, irreducible, Birkhoff integrable Hamiltonian 
system is isomorphic to one of the following diagrams: 

o 1 o 1 




^1 1 1 

o — o o — o 



6 



6 



(a) 



(b) 



o 

(d) 



(e) 



1 1 1 2 2 

o — o — — o 



2 2 2 1 

o — o — o=o 



(f) 



(g) 





Remark 1 In the list of diagrams we have omitted some cases that occur 
as sub-graphs of diagrams (a)-(k) (by truncating one or more vertices). 
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Remark 2 The Dynkin type diagram determines only the angles be- 
tween pairs of vectors in A. In order to reconstruct the ratios of lengths 
of vectors in A we assign to the zth vertex a coefficient proportional to 
the square of the length of Vj. This explains the numbers appearing on 
the vertices of the diagrams. 

We have to stress that this classification gives only necessary con- 
ditions for a system of type (j2J) to be Birkhoff integrable. The inte- 
grability for each system in the list should be established case by case. 
As we already mentioned, the integrability of systems (a)-(g) was es- 
tablished in [1]. The solution of these generalized periodic Toda 
lattices (associated with affine Lie algebras) was obtained by Good- 
man and Wallach in [T3]. The graph (i) corresponds to a Hamiltonian 
system in two degrees of freedom with potential 

e qi + e q2 + e- qi ~ g2 + e"(^) . 

The additional integral can be found in |16j . 

Sklyanin [23] pointed out another integrable generalization of the 
Toda lattice: 
(9) 

n j n—1 

H{q u . . . ,q n , ft, . . . ,p n ) = 2^+1] e ft - 9i + 1 +a 1 e 9l +/V 91 +ane~ <z ™+/3ne- 

i=i i=i 

He obtained this system by means of the quantum inverse scattering 
R-matrix method. This corresponds to system (h) with associated 
Hamiltonian The case n = 2 corresponds to the potential 

V = e qi ~ q2 + Cl e 2q2 + c 2 e q2 + c~ 91 + c 4 e" 2gi . 

Annamalai and Tamizhmani [3] demonstrated the integrability of this 
particular case by using Noether's theorem. The second integral is of 
fourth degree in the momenta. A special case of was considered in 
[8] using the generalization of KM (Volterra) system due to Bogoyavlen- 
sky. The case n = 3 (as well as the general case) is treated in Ranada 
|22j . Ranada proved integrability by using a Lax pair approach. The 
additional integrals are of degree 4 and 6. Ranada's approach will be 
used in proving the integrability of system (CQ). It is generally believed 
that system (k) in non-integrable. In fact the arguments in [16] support 
the non- integrability of system (k). This leaves only system (j) which 
is treated in the present paper. 

The only essential progress to our knowledge in proving the integra- 
bility of system (CQ) are the two papers [9] and pH] • The case n = 4 was 
studied in [TH] . Based on calculations of the Kovalevskaya exponents it 
was conjectured that the degrees of homogeneity of the additional first 
integrals are 4, 6, 8. Our results confirm that prediction. We reproduce 
here a small part of the table from [10] which contains the Kovalevskaya 
exponents corresponding to the various values of the indicial locus c: 
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Vector c 


KE exponents 


(0,6,0,10,6,0) 
(3,4,3,3,0,0) 
(7,12,0,15,0,16) 
(8,14,0,18,10,0) 
(1,0,1,4,3,0) 
(3,4,3,0,1,0) 
(7,14,0,18,10,0) 


-5,-3,-2,-1, 1,1,1,2,4,6,8,8 

-3,-2,-1, 1, 1,1,2,3,4,5,6,8 

-7,-5,-3,-1, 1,1,1,2,4,6,8,14 

-7,-5,-3,-2,-1,1,1,2,4,6,8,16 

-3,-2,-1,-1,-1,1,1,2,2,2,4,8 

-3,-2,-2,-1,-1,1,1,2,2,3,4,8 

-7,-5,-3,-2,-1,1,1,2,4,6,8,16 



In a brief communication by Emelyanov [H] it was reported (without 
proof) that the system ([T]) was integrable for n = 4. In fact a few 
leading terms of the integrals were displayed but our leading terms of 
the integrals do not match exactly those results. However, it is possi- 
ble that the integrals calculated by the author of [9] are combinations 
of the integrals of the present paper. To quote from [9]: "A study 
of Kovalewski indices shows that it is possible to assume for n = 5 
the degrees of additional integrals are 4, 6, 8, 10 etc. If integrability is 
established for an arbitrary n, then one can state that only systems 
whose Dynkin schemes are isomorphic to those of Fig. 1 except for 
scheme (k) are Birkhoff integrable" . 

2. D n TODA SYSTEMS 

Since the Lax pair of system (CQ) uses in an essential way the Lax pair 
of the D n Toda lattice we include a review of this Bogoyavlensky-Toda 
system following [5], [6], [7j. 

The Hamiltonian for the D n Toda lattice is 

1 n 

(10) H = - 22 pj + e qi ~ q2 + ■■■ + e qn - 1 - qn + e qn ~ 1+qn n > 4 . 

i 

We make a Flaschka-type transformation, F : R 2n — > R 2n defined 

by 

F : (gi, ■ ■ -,q n ,Pi, ■ ■ ■ ,Pn) -> i a i, ■ ■ ■ ,a>n,W, ■ ■ ■ A) , 

with 

(11) oj = I e ^ ( « +l) , i = 1,2,. a n = - e ^ qn - 1+qn) , 

bi = --pi, i = l,2,...,n. 

Then 
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(12) 



a„ 



cti (b i+1 - hi) 
-a n (b n -i + b n ) 
2 (a 2 - aU) 
2{a 2 n + a 2 n _ 1 1 



1,2,. ..,n 



1, 2, . . . , n — 2 and i 



n 



n-2) 



These equations can be written as a Lax pair L = [B,L], where L is 
the symmetric matrix 



(b\ ai 
a\ ' • ■ 



\ 





O'n-1 


-a n 





O n -l 


K 





O n 


—On 





~b n 


—o n -i 







— O'n-1 





(13) 



—0\ 

\ -ax -b x J 

and B is the skew-symmetric part of L (In the decomposition, lower 
Borel plus skew-symmetric). 

The mapping F : R 2n — > R 2n , (q^Pi) -> (o^fy), defined by (ITTD . 
transforms the standard symplectic bracket into another symplectic 
bracket ttx given (up to a constant multiple) by 





{ai,bi} = 


o Oi 


2=1,2,. 


. ,n 


(14) 


{ai,b i+1 } = 


ia 


2=1,2,.. 


. , n — 1 




{o n ,b n -i} = 


z i 

2 • 







We obtain a hierarchy of invariant polynomials, which we denote by 

H2, H4, . . . , H2n, ■ • • 

f Tr L 2i . 

polynomials are 2, 4, . . . , 2n — 2. We also define 



defined by i^i = 27 Tr L 2 \ The degrees of the first n — 1 (independent) 



P n = VdetL . 

The degree of P n is n. The set {H 2 , i7 4 , . . . , H 2n _2,P n } corresponds 
to the Chevalley invariants for a Lie group of type D n . 
Taking 



Ho 



1=1 



i=l 



as the Hamiltonian we have that 

7TiV# 2 

gives precisely equations (|T2l . 
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3. QUADRATIC LAX PAIRS 



We employ the method used by Ranada in [22]. The key idea is the 
following: If (L, B) is a Lax pair, so is (L 2 , B). This follows easily from 



d 
dt 



L 2 = [B,L]L+L[B,L] = (BL—LB)L+L(BL—LB) = BL 2 —L 2 B = [B,L 2 
In the case of Da Toda lattice, L 2 is the matrix 



(15) 

/ al + bl 
ai(6i + 6 2 ) 









V 



ai(6i + 6 2 ) 
a{ + al + bl 
a 2 (6 2 + 63) 
02 03 
—0,20,4 






2 , 2 
a 2 + a 3 



0\02 

02(62 + 63) 

al + bl 
03(63 + 64) 
a 4 (6 4 - 63) 
20304 







0203 
03 (63 + 64) 
a'l + al + b'i 

— 20304 
a 4 (&4 - 63) 
—0204 






—0204 
a 4 (&4 - 63) 
—2a 3 a4 
al + bl 
03(63 + 64) 
0203 




al 



and 



(16) 



B 



( 


ai 

















° \ 


—Ol 





02 




















—02 





03 


— 04 

















-a 3 








a 4 














a 4 








-a 3 

















— 04 


03 





-02 

























-01 


V 

















a 1 


/ 






20304 
04(64 — 63) 
03(63 + 64) 
al + a 2 3 + al + b 2 3 
02(62 + 63) 

Ol02 



Note that the Lax equations L 2 = [B, L 2 ] gives the same equations 
( IT2|) plus some consistency conditions. 

4. The case n = 4. 

We now consider the special case of ([T]) for n — 4. The general case 
follows easily from this special case. Let 

4 3 

(17) # = J2 2 + e9i_<?i+1 + e<?3+94 + e ~ 91 + 

i=l i=l 

We make a Flaschka-type transformation, F : R 8 — 



-2qi 







—0204 
0203 
02(62 + 63) 

of + ol + &| 
Ol(6i + 6 2 ) 



\ 









Ol02 

01(61 + 6 2 ) 
o? + 6f / 



R 9 defined by 



F : (q 1 ,...,qA,p ll ...,p i ) 



with 
(18) 



e 2 



hili-H+i) 



ai, 



,a 5 ,b 



i- 



1 

— < 

2 



i = 1,2, 3, 04 = -e 2 



1 



5(93+94) 



k = -^Pi, 1 = 1,2,3,4. 

We obtain the following equations of motion: 



64) 



«5 



v/2 



e 2 



91 
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ai(b i+1 -bi) i = 1, 2, 3 

-a 4 (6 3 + h) 

2a? -a 2 - 4a* 

2 (a| — a?_ x ) i = 2 and 2 = 4 

2(a^ + a| — a|) . 

Note that the Hamiltonian in the new variables takes the form 

4 4 

(20) F = J]6? + 2^g 1 2 + ^ + 2^. 

i=l i=l 

The image of the symplectic bracket is now the following extension 
of bracket (fT4"|) in the new phase space (ai, . . . , 05, b±, . . . , 64). 

{a-iA} = -iai i = 1,2,..., 4 

(21) { a5 ' bl \ = ! as . 1 . . 

{a h b i+1 } = i = 1,2,..., 3 

{04,63} = — -^a^. 

We denote this bracket by W\. Using the Hamiltonian (1201) and the 
above bracket W\ gives equations (TT9"j) as is easily checked. 

We now present a Lax pair for equations (|T9l) . The method of ob- 
taining this Lax pair is based in the ideas of [22] and we omit the 
details. 

The Lax pair has the form A = [C, A] where the matrices C and A 
are perturbations of the matrices B and L 2 respectively. More precisely, 
A is the same as L 2 except that 

A u = a 2 + b\ + a\ + 2a\ = A 88 , 
A 12 = a 1 (b 1 + b 2 + V2iaj) = A 87 , i 2 = -1 

A 2 i = ai(&i + b 2 - V2ial) = A 78 . 

On the other hand C differs from B only at two diagonal positions. 
It is the same as B except that 

Cu = V2ial = C 88 . 

It is a simple calculation to show that the equations A = [C, A] are 
a matrix form of equations ( |T9l) . 
Define h 2i = \ Tr A i i = 1, 2, 3, 4. 

1 4 4 

fca = -Tr A = H = ^ b 2 + 2 ^ a 2 + a* + 2a* . 

i=i i=i 



(19) 



d 4 
as 

&, 
6 3 



Let 
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ti = a\ + \ai + at 



1 ~ 2 5 ~ 5 

t 2 = a\ + a\ 

t 3 = a 2 2 + a\ + a\ 

^4 — flg -|- O4 . 

Then /i 4 = |Tr ^4 2 is given by 

4 

/u = 6j + ^ + ^ + 6i + E 4t ^ 2+ 

1=1 

4a 2 6i6 2 + 4:dlb 2 b 3 + A(a 2 3 - 04)6364 + s(a 1 , a 2 , a 3 , a 4 , a 5 ) , 
where 

s = Aa\al+Aalal+Aalal+l2alal+M\al+2^ 

Similarly h 6 has the form J2t=i + J2t=i 12£»&f+ other terms even 
in the momenta. 

Finally fog has the form £) i=1 6^ + ^2 i=1 8tj6f + other terms even in 
the momenta. 

The set of functions {/12, h^, he, h$} provide a maximal set of inde- 
pendent integrals in involution. The fact that the leading terms of the 
integrals are functions of the momenta only, agrees with the general 
facts proved in [16J. 

5. The general case 

The procedure for the general case is very similar to the case n — 4. 
There is no interaction between q\ and the newly introduced variables, 
therefore the procedure is quite identical. 

Consider the Hamiltonian 

n 1 n— 1 

H = -Pi + ^2 e " +1 + e qn - 1+q " + e~ qi + e~ 2qi . 

i=l i=l 

We make a Flaschka-type transformation, F : R n — ► R n+1 defined 

by 

F : (q l ,...,q n ,p 1 ,...,p n ) -> (01, . . . , a n+1 , h, . . . , b n ) , 

with 



(22) 



^ Mli-li + l) A — 1 O n _ l„^(9n-l+gn) ' 



1 



ai = -e^«-'^, z = l,2,...,n-l, a n = -e^-™ , a n+1 = -=e~2 

6< = ~\pi-> i = 1,2, . . .n . 

We obtain the following equations of motion: 
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dj = cii(b i+1 -bi) i = 1,2, ...,n- 1 

«n = -a n {b n -i + b n ) 

, OQ x <Wl = Un+lbl 

' ' 6i - 2r7 2 -o 2 -An 2 

6j =2 (a 2 — a 2 _ x ) z = 2, 3, . . . , n — 2 and i = n 

b n -i = 2(a 2 + a 2 _! - a 2 _ 2 ) . 

Note that the Hamiltonian in the new variables takes the form 

n n 

(24) H = Y,b* + 2j2a-+al + i + K + i- 

i=i i=i 

The image of the symplectic bracket is now the following extension 
of bracket ( TUE|) in the new phase space (ai, . . . , a n+ i, 6i, . . . , 6 n ). 



HA} = 




i = l,2,.. 


. ,n 


{a n+ i,6i} = 


9 u ! 






{di,b i+1 } = 


2=1,2,.. 


. , n — 1 


{a n ,6 n _i} = 


^ 1 

2 ®n ■ 







We denote this bracket by W\. Using the Hamiltonian f[2Tl) and the 
above bracket w\ gives equations ( 1231) as is easily checked. 
We now obtain a Lax pair for equations ( 1231) . 

The Lax pair has the form A = [C, A] where the matrices C and A 
are perturbations of the matrices B and L 2 respectively. More precisely, 
A is the same as L 2 except that 

A u = al + b\ + a 2 n+1 + 2a* +1 = A nn , 
A 12 = ai(6i + b 2 + V2ial +1 ) = A {n+l)n , i 2 = -1, 

A 21 = ax{bx + b 2 - V2ia 2 n+l ) = A n{n+1) . 

On the other hand C differs from B only at two diagonal positions. 
It is the same as B except that 

C\\ = v^id^^ = C nn . 

It is a simple calculation to show that the equations A = [C, A] are 
a matrix form of equations (1231) . 

Define h 2i = § Tr A i % = 1, 2, . . . , n. 
e.g., 



h 2 = ^Tr A = H = jS 2 + 2^a 2 + a 2 n+1 + 2a 4 n+1 , 

i=l i=l 



2 

and h 2 i is a homogenous polynomial in the bi of degree 2i. 
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